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Abstract We discuss subtle points of the momentum balance for radiating particles 
in flat and curved space-time. An instantaneous balance is obscured by the presence 
of the Schott term which is a finite part of the bound field momentum. To establish 
the balance one has to take into account the initial and final conditions for acceler- 
ation, or to apply averaging. In curved space-time an additional contribution arises 
from the tidal deformation of the bound field. This force is shown to be the finite 
remnant from the mass renormalization and it is different both form the radiation re- 
coil force and the Schott force. For radiation of non-gravitational nature from point 
particles in curved space-time the reaction force can be computed substituting the 
retarded field directly to the equations of motion. Similar procedure is applicable to 
gravitational radiation in vacuum space-time, but fails in the non-vacuum case. The 
existence of the gravitational quasilocal reaction force in this general case seems im- 
plausible, though it still exists in the non-relativistic approximation. We also explain 
the putative antidamping effect for gravitational radiation under non-geodesic mo- 
tion and derive the non-relativistic gravitational quadrupole Schott term. Radiation 
reaction in curved space of dimension other than four is also discussed 



1 Introduction 

One of the major tasks of gravitational wave astronomy is the precise theoretical pre- 
diction and observational measurement of gravitational waveforms from an inspiral 
fall of compact bodies into the supermassive black hole. This requires knowledge of 
the orbits with account for gravitational radiation reaction. Radiation gives rise to 
the reaction force which can be incorporated into the equations of motion. The stan- 
dard strategy to get the reaction force consists in substitution of the retarded field 
produced by the body into its equation of motion. The resulting equation is believed 
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to give a correct description of the instantaneous effect of radiation on the motion 
of the body. If the effect of reaction is small with respect to the main force, one 
can treat it adiabatically. But when the reaction force is not small, the situation is 
more subtle: the balance equations involve not only the kinetic energy-momentum 
of the body and that of radiation, but also a variable contribution of the bound field. 
Generically, an instantaneous loss of the energy-momentum by the body is not equal 
to the energy-momentum carried away by radiation. 

This situation looks particularly simple in the flat space Maxwell-Lorentz elec- 
trodynamics, where due to linearity of the equations one can use the notion of the 
point-like particle described by the delta-function. The Lorentz-Dirac equation ob- 
tained by substituting the retarded field into the equation of motion and performing 
the mass renormalization is expected to describe the motion of the particle with 
account for radiation loss. However, it turns out that the momentum loss due to ra- 
diation gives only a part of the reaction force. Although the initial system of the 
charges and the Maxwell field obey the overall momentum conservation equations, 
the Lorentz-Dirac equation, treated as the particle equation of motion, violates the 
naively expected balance between the particle momentum and the momentum of ra- 
diation. The difference is given by the so-called Schott term, which is the third order 
total derivative of the coordinate. A careful comparison of the Dirac derivation 0~1|2) 
and the Rohrlich analysis J3]|4][5) nas l e d Teitelboim (6) (see also Q) to interpret 
the Schott term as the finite contribution of the bound (non-radiated) momentum of 
the charge remaining after the mass renormalization (an explicit proof of this claim 
was given in [8 1). This may look strange, since we are used to think about the bound 
field as a stable Coulomb coat which is spherically symmetric in the instantaneous 
rest frame or pancake shaped for a relativistic charge. But this simple picture is valid 
only for constant velocity. Once acceleration is non-zero, the energy-momentum 
carried by the Coulomb coat becomes variable. The most surprising fact is that this 
momentum is simply proportional to the acceleration, and thus its derivative (the 
force) is the third derivative of the particle coordinate. Of course, the split of the 
total field into radiation and the bound field has to be done at any distance from the 
charge, not only in the wave zone. This was given by Rohrlich \5], for more recent 
discussion see the book by B. Kosyakov |9|. It is worth noting that in this case the 
finite part of the reaction force is entirely given by the time-antisymmetric part of 
the particle field (half-difference of the retarded and advanced potentials). 

Therefore the energy-momentum balance of the system consisting of the accel- 
erated charge and its Maxwell field includes three, but not just two, ingredients: the 
particle momentum, the momentum carried by radiation, and the bound electromag- 
netic momentum. The radiation momentum can be extracted both from the particle 
momentum and indirectly from the bound momentum. This explains the origin of 
radiation of the uniformly accelerated charge, in which case the total reaction force 
is zero and thus the kinetic particle momentum is constant. While the charge is un- 
dergoing a constant acceleration, its bound electromagnetic momentum decreases 
and is transferred to radiation. Physically, however, the acceleration has to start at 
some moment and to finish at some moment, and during the stages of acquiring and 
loosing the acceleration the bound momentum is exchanged with the kinetic mo- 
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mentum. Therefore, the total energy-momentum loss of the charge will be equal to 
the momentum carried away by radiation. But an instantaneous balance is obscured 
by the presence of the Schott term. Another simple situation is periodic motion. 
Since the ambiguous Schott term is total derivative, its contribution vanishes if one 
integrates over the period, or, equivalently, averages over the period. 

For radiation of linear fields of non-gravitational nature or of the linearized grav- 
itational field in curved space-time the situation is more complicated because there 
are no local conservation laws unless the Killing symmetries are present, and be- 
cause the tail term can not be found in the closed form ifTUl [TP (for a review see 
fl2l . Another new feature in this case is that the reaction force contains a finite time- 
symmetric contribution coming form the half-sum of the retarded and advanced 
potentials - the force found by DeWitt-DeWitt |[T3l in linearized gravity and later 
rediscovered in the full General Relativity by Smith and Will HH (the WWSM 
force). This force, however, has nothing to do with radiation and so the work done 
by this force is not expected to contribute into the energy balance for radiation. In 
the stationary space-time the total energy of the particle and the field is conserved, 
so one can expect that the energy balance between radiation and the particle en- 
ergy loss will hold integrally or in average. In the case of axial symmetry, similar 
considerations apply to an associated angular momentum. An explicit proof of this 
balance for scalar, electromagnetic and linearized gravitational radiation in the Kerr 
space-time was given in 03] ■ Apparently this work was not properly understood 
and was criticized in a number of papers until 2005, when analogous calculations 
were performed by other people (for a review and further references see lTT6l ). 

The case of linearized gravity is not entirely similar to other linear fields in the 
curved background, however. In fact, linearized gravity on the non-vacuum back- 
ground is not a consistent linear theory since the full Bianchi identities do not allow 
for the harmonic gauge necessary to locally disentangle the linearized Einstein equa- 
tions fm . Physically this means that the proper source of the gravitational radiation 
from the point particle is not just its energy-momentum tensor, but it also includes 
contribution from the perturbed source of the background. Thus the total source is 
non-local, which raises doubts that there might exist a local equation describing ra- 
diation reaction (here we mean the non-locality stronger than that of the tail term, 
which is still localized on the world-line of the point particle). 



2 Energy-momentum balance equation 

Consider the interacting system of N point charges and the Maxwell field in 
Minkowski space-time which is described by the coupled system of equations 




£1=1 




(1) 
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It has 3N + °° degrees of freedom, where °° stands for the Maxwell field. The corre- 
sponding energy-momentum conservation equation is 

d v {T flv +T^)=0, (2) 

where 

(3) 

We have introduced the bare masses in anticipation of mass-renormalization. 
Since the Maxwell equation is linear, one can decompose the total field in the vicin- 
ity of any given charge e a into the sum of the field generated by the other N —1 
charges, F^ t (regular at its location) and the retarded field of e a , F^ eV In spite 
of the fact that the total field acting on e a , Fa V = F^ + F^ et diverges at x^ = Za, 
energy-momentum conservation is ensured by the equations of motion. The required 
mass renormalization does not change this statement. 

The field F^ et describes radiation and the bound field, which both contribute to 
the field energy-momentum. The overall conservation equation does not distinguish 
between these two parts, so an additional analysis is needed. The retarded potential 
at a given point x of space-time depends on the world-line variables taken at the 
moment of proper time s ret (x) defined as the solution to the equation 

#%=0, fl^^-Z^e,), (4) 

satisfying x° > z°- The advanced solution to the same equation with z° > x° refers 
to the advanced proper time s a i v (x). Introducing the invariant distance 

p=v M (* ret )^, v" = ^, (5) 

which is equal to the spatial distance |R| = |x — z(5 , ret )| between the points of emis- 
sion and observation in the momentarily co-moving Lorenz frame at the time mo- 
ment x° = z°{s let ), one can present the retarded potential as (we omit the index a): 



^ret( x ) — 

r 



(6) 

SretW 



Introduce the normalized null vector c^ 1 — R^/p, such that vc = 1, and the unit 
space-like vector = — , u 2 = — 1 . The following differentiation rules then 
hold: 

c M =d ll s m (x) 7 <? M p = v M +Ac M , <9 M c v = -(5^-v M c v -c M v v -Ac M c v ,) (7) 

P 



where A = p = p (ac) — 1 . The retarded field strength will read: 
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F ,v = e (p(*c)-l) v[ e { (g) 
P P 

The retarded potential in Minkowski space admits a natural decomposition with 
respect to T-parity: 

<t=<if+<d: (9) 

where the radiative part A^ ad = j (A^ t — A^ dv ) obeys an homogeneous wave equa- 
tion, while the self part A^ elf = \ (A^, t +A^ dv ) has a source at x = z(s). One could 
expect that only T-symmetric A^ elf corresponds to the bound field, but it is not so. 
For an accelerated charge the situation is more subtle. 



2.1 Decomposition of the stress tensor 

Constructing the energy-momentum tensor T^ v with the retarded field F r ^ t , one 
finds that it admits a natural decomposition: 

1 - 1 cmit+ 1 bound' 



where the first term is selected by its dependence on p as p 



-2. 



£ M v . ((ac) 2 + a 2 )c^c v 

* emit n 2 ' ^ ^ 



while the second contains higher powers of p 



-i. 



F 



_ a^c v > +2(ac)c^c v - (ac)v^c v > v^c v ) - c^c v - t/ v ^ 

^bound — ^3 ' ^4 ' ' 

where symmetrization without 1 /2 is understood. 
The "emit" part (fTTT i has the following properties: 

• it is the tensor product of two null vectors c^, 

• it is traceless, 

• it falls down as |x|~ 2 when |x| — > °o, 

• as follows from the differentiation rules (|7}, it is divergence-free without assum- 
ing the validity of the equations of motion: 

*A = 0. (12) 

All these features indicate that describes the outgoing radiation. 

Since the total energy-momentum tensor including the contribution of charges is 
(on shell) divergence free, with account for ( fT2l we find that 



6 



Dmitri Gal'tsov 



^ft o v und + ^f>=0, (13) 

so the bound field momentum can be exchanged with the particle momentum. Note, 
that outside the world-line the bound stress tensor is divergence-free. It is also worth 
noting, that Eq. (fT2l i does not mean that there is no reaction force acting on a particle 
which counterbalances the emitted momentum. 

Consider now the total balance of forces. The conservation of the total four- 
momentum (O implies that the sum of the mechanical momentum and the momen- 
tum carried by the electromagnetic field is constant (for simplicity we do not include 
the external field): 

jgmgch ! d Pem = Q (U) 
ds ds 

Here the mechanical part is proportional to the bare mass of the charge 

<ech = /^v=mV, (15) 
while the field part is given by 

pL = j fv v dZ v , (16) 

where integration of the electromagnetic stress tensor is performed over a space-like 
hypersurface whose choice will be specified later on. It has to be emphasized that 
the stress tensor of the electromagnetic field is constructed in terms of the physical 
retarded field. According to the above splitting, we can write 



mech 



/emit "'"/bound' (XT) 



/emit - ~ / ^emit^V, /hound ~ /^hound^V- (18) 



ds 

emit^ v ' 

On the other hand, the derivative of the bare mechanical momentum of the charge 
can be found substituting the retarded field in the equation of motion. In this case it 
is useful to decompose the retarded field according to (0, obtaining another split of 
the mechanical momentum: 



Z£mech_ pMV _.( F HV p nv\ _ M , M (]Q s 



^mech _ „^iv„ _ „ rcMV, ipliV\ _M . M 



Now, somewhat unexpectedly, / r ^ d ^ fj^ t and f^ f ^ f£ ound , the difference being 
called the Schott term 0: 

/ rad = / smit + A Schott ' /elf ~ /bound — / Schott " 

Clearly, 

/elf + /rad ~ /bound + /emit ' (21) 
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as expected. Note that both f£ ound and contain divergences which mutually can- 
cel in Eq. ( f2Tb . 

The forces /^ lf and f^ d can be found using the Green functions [2] 

G self (Z) = 5(Z 2 ), G rad (Z) = |J5|5(Z 2 ), (22) 

where ~Z tl {s,s') =z^(s) — z'^(s'). Substituting the value of the electromagnetic 
field generated by the charge on its world-line one obtains 

f{s)=2e 2 J Z^{s 1 s')v v \ S ')v v {s)^ i G{Z)ds' ) (23) 

for both and Due to delta-functions, only a finite number of Taylor expan- 
sion terms in a = s — s 1 contribute to the integral. In the four-dimensional case it is 
sufficient to retain the terms up to a 3 : 

2Z^(s,s')v v \s')v v {s) =v^a 2 -^+vV)(7 3 + <9(cT 4 ). (24) 

Taking into account that Z 2 = a 2 + 0(a 4 ), the leading terms in the expansions of 
derivatives of the Green functions will be 

&g*M = ±*&), ^ G ^ = Mw\ 5[a2) )- (25) 

Regularizing the delta functions of a 2 by point- splitting 

8(a )= hm S(a -e )= hm — — — -. (26) 



fH* = -^ /^d = ^-(v^ 2 + ^). (27) 



with a prescription that the limit should be taken after evaluating the integrals, one 
finds 

2E ' /rad_ 3 

After the mass renormalization, mo + = m we get the Lorentz-Dirac equation 

ma> i = 7 ^-{v^a 2 + d^). (28) 

The first terms at the right hand side is equal to the derivative of the momentum 
carried away be radiation, 

/emlt ds 3 K ' 
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Its independent evaluation by integration of T^L, can be found in (8). The second 
total derivative term is the Schott term. It is worth noting, that within the local cal- 
culation, the Schott term originates form the T-odd part of the retarded field. 



2.2 Bound momentum 

An explicit evaluation of the bound momentum associated with a given moment of 
proper time s on the particle world-line, 

PL»= / fc,nd^v (30) 

was given in J8], which we follow here. First of all one has to choose the space-like 
hypersurface E(s) intersecting the world line at x^ =7^{s). A convenient choice 
will be the hyperplane orthogonal to the world-line 

v m (j)(j^-^(j))=0. (31) 

The integral ( |30l is divergent on the world line. We introduce the small length 
parameter e, the radius of the 2-sphere dY £ (s) (Fig. 1), defined as the intersection 
of the hyperplane (|3~T1 i with the hyperboloid (x — z{s)) 2 = — e 2 . We also introduce 
the sphere 3Yr(s) of the large radius R defined as the intersection of E(s) with 
the hyperboloid (x — z{s)) 2 = —R 2 . The total field momentum can be obtained as 
the limit e — >• 0, R — > °° of the integral over the domain Y(s) C E(s) between the 
boundaries dY £ (s) and 3Yr(s). 

Let us evaluate the variation of this quantity between the moments s\ and S2 of 
the proper time on the world-line of the charge 

ApL= J T^ v dE v - J T^ v dE v . (32) 

Y( S2 ) r(n) 

For the bound momentum it is convenient to consider the tubes S £ and Sr formed 
as sequences of the spheres dY £ (s) and 3Yr (s) on the interval sG sz] and to 
transform this quantity to 

^Pbound = / T tZxd dS V - J T bZnd dS V ( 33 > 

Sr S e 

in view of the conservation equation for T^ nd outside the world-line (see the re- 
mark after Eq. (fT3Tl). Here the integration elements dS v are directed outwards to the 
world-line. The contribution from the distant surface Sr vanishes if one assumes that 
the acceleration is zero in the limit s —> — °° [6|. This is non-trivial: though the stress 
tensor (fT2l decays as R~ 3 at spatial infinity, the corresponding flux does not vanish 
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Fig. 1 Integration of the bound electromagnetic momentum. Here E(si) is the space-like hyper- 
plane transverse to the world-line z^(s) intersecting it at the proper time S\ (similarly Ufa)). The 
hypersurfaces S £ and Sr are small and large tubes around the world-line formed by sequences of 
the 2-spheres dY e (s) and 9Yr(s) for s 6 [si, st\. The domain Y(s%) C (similarly Y(s\ )) is the 

3-annulus between 9Yr(s2) and dY e (s2). 



a priori, because the surface element contains a term (proportional to the accelera- 
tion) which asymptotically grows as E? . As a consequence, the surviving term will 
be proportional to the acceleration taken at the moment s mi of the proper time, where 
■Sret — > —°° in the limit R — > °o. Finally we are left with the integral over the inner 
boundary S £ only. To find an integration measure on S £ we foliate the space-time do- 
main shown in Fig. 1 by the hypersurfaces E (s) parameterized by the spherical coor- 
dinates r, 0j = 0, 02 = <P- Introducing the unit space-like vector n^(s, 0,), n^n^ = — 1 
transverse to v^, we use the coordinate transformation = (s) + rn >1 (s, 0,). The 
induced metric on S £ reads dS^ — e 2 [l — e(an)]n^dsd£2, and hence 

Abound f - F 



ds 



[ e 2 [l-e(an)]T^ m<i n v dn, (34) 
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where the limit £^0 has to be taken. One has to expand T^ ound in terms of e. 
In fact, the energy-momentum tensor depends on the space-time point x^ through 
the quantity p, depending directly on x 11 , and also through the retarded proper time 
■Sret- We have to express the resulting quantity as a function of the proper time s 
corresponding to the intersection of the world-line with the space-like hypersurface. 
We write T^° and in terms of the null vector R^ = c^p: 



An MV _ a { ^R v) (2(aR)-l)RVR v (1 - ( a R)) y^R v ) ^ 

Abound - —^T + ~6 + p5 ^4 ■ 



and expand as 

R H =xi x _ z M( iret ) = en /i + v H a _ 1 M CT 2 + l ^ a i + ff f^ (36) 

2 6 

where a = s — s ra > and all the vectors are taken at s. This expansion in powers 
of o has to be rewritten in terms of £. The relation between the two can be found 
from the condition R~ = 0: 

3 

a = e + y e 2 + (9(an) 2 + a 2 - 4dn) + 6\e A ). (37) 
Substituting this into Eq. (l36i l and further to d35l l we find: 

2 ^ 2 

^ - ^ I** { + + [(M 2 + fl2 /3) 

+ ({an) 2 +a 2 /2)n >i -2d fi /3 + 3(an)a^/4]}d£2. (38) 

The leading divergent term 1 /e 2 disappears after angular integration. Thus we ob- 
tain i): 



f H _ Abound _ ^ . 2e 2 M 
/bound- ^ - 2 £ 3 



Here the first divergent term has to be absorbed by the renormalization of mass, 
while the second is the finite Schott term. Comparing this with ( |27| | we confirm the 
identity ( f2Tb . Note that a priori the regularization parameter e here (the radius of 
the small tube) is not related to the splitting parameter of the delta-function in the 
previous local force calculation. But actually they give the same form for divergent 
terms, for which reason we use the same symbol e for both of them. With this con- 
vention, the divergent terms in the momentum conservation identity (f2Tb mutually 
cancel. 

The significance of the Schott term in the balance of momentum between the ra- 
diating charge and the emitted radiation is not always recognized in the text-books 
on classical electrodynamics. Instead, its presence is often interpreted as a drawback 
of classical theory, since formally it may lead to self-accelerating solutions. Mean- 
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while such solutions must be discriminated as unphysical since they do not satisfy 
proper initial/final conditions which should be imposed on the third order equation 
of motion |2|. From the above analysis it is clear that the Lorentz-Dirac equation, 
although formulated in terms of particle variables, actually describe the composite 
system consisting of a charge and its bound electromagnetic momentum. The redefi- 
nition of the particle momentum joining to it the bound electromagnetic momentum 
obscures the problem of interpretation. It is better think of the Schott term as the 
field degree of freedom and interpret the Lorentz-Dirac equation as the momentum 
balance equation for the total system including the electromagnetic field. 

An instantaneous momentum balance is not just the balance between the particle 
and radiation, the energy-momentum can be also transferred between the particle 
and the field coat bound to it. Or, the radiated momentum is not always taken from 
the mechanical momentum of the charge, but by virtue of the Schott term, it can be 
extracted from the field coat too. This is what happens in the case of the uniformly 
accelerated charge, when the total reaction force instantaneously is zero, while ra- 
diation carries the momentum away at a constant rate. The balance in ensured by 
the Schott term. However, the constant acceleration can not last infinitely long. One 
has to consider the switching on/off processes in order to understand that finally the 
energy-momentum of radiation is taken from the particle. This consideration clari- 
fies the necessity of the time averaging or integration over time needed to establish 
the momentum balance between radiation and the source particle. The equation of 
motion including the reaction force instantaneously does not imply the equality of 
the radiative momentum loss and the particle momentum. This feature is general 
enough, it is also applicable to radiation of non-gravitational nature from particles 
moving along the geodesies in curved space-time, as well as to the gravitational 
radiation. 



2.3 The rest frame (non-relativistic limit) 

In the rest frame of a charge the recoil force has no spatial component. This is due 
to the fact that radiation in two opposite directions is the same so that the spatial 
momentum is not lost by radiation, though the energy is lost. Hence the total spatial 
component of the reaction force is presented by the Schott term, namely 

fschott - ^ 2 a. (40) 

The work done by this force 



/ 



/ ^e 2 a • vdt — — J ^e 2 a 2 dt + boundary terms (41) 



correctly reproduces the radiative loss in the rest frame. (Boundary terms should 
vanish by appropriate asymptotic switching on/off or periodicity conditions). 
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3 Flat dimensions other than four 

Recent interest to models with large extra dimensions motivates the study of ra- 
diation and radiation reaction in dimensions other than four. It turns out that the 
radiation picture is substantially different in even and odd dimensions because of 
the different structure of the retarded Green's functions for massless fields in the 
coordinate representation 1 18 1 (they still look similarly in all dimensions in the mo- 
mentum representation). In even dimensions the retarded potential is localized on 
the past light cone (Huygens principle) so the situation is qualitatively similar to 
that in the 4D case. In odd dimensions it is non-zero also inside the past null cone 
though radiation still propagates along the null rays. In 3D, for instance, the scalar 
Green's function reads 

G%(X) = &(X°)&(X 2 )(X 2 )(- l / 2 \ X^=x^-x^. (42) 

It does not-contain the "direct" part singular on the light cone. Green's functions in 
higher odd dimensions D = 2n + 1 can be obtained by the recurrent relation |[T8l[T9l 

j pin- 1 

G^ + \X)~^- (43) 

In particular, in 5D 

both the direct and the tail parts are present. It turns out that the direct part regular- 
izes the tail contribution to the field stress proportional to the derivative of G which 
otherwise would be singular outside the world line. 

In even dimensional space-times the split of the retarded potentials into the time 
symmetric and the radiative parts leads to purely divergent self-force and a finite 
radiative part: 

Ci = fL C= fini te. (45) 

Divergent terms are Lagrangian type and can be absorbed by introducing suitable 
counterterms. Since the Coulomb dependence is more singular at the location of 
the source in higher dimensions, the self-action gives rise to larger number of di- 
vergent terms 1/e", where n is changing from unity to the integer part of D/2 — 1. 
The highest divergence is absorbed by the renormalization of mass, while to absorb 
other divergencies the counterterms are needed depending on higher derivative of 
the velocity. These are not present in the initial action, so higher-dimensional clas- 
sical theories are not renormalizable. In 6D one has two divergent terms (which in 
terms of the field split correspond to / se if fZUl : 

1 1 /3 3 1 \ 

fL = - ^3 a v + Ye [ 4 V/i (da) + 8 fl2flA1 + 4 &VL ) ' (46) 
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the leading being eliminated by the mass renormalization and the subleading requir- 
ing the counterterm |j20l 

Si = -^ 1} Jdfds, (47) 

which leads to the Frenet-Serret dynamics lETI unless the renormalized value K^' = 
0. For each two space-time dimensions one new higher-derivative counterterm is 
needed to absorb divergencies. 

The split of the field stress-tensor built with the retarded field into the sum of the 
emitted and bound terms is also possible in all even dimensions, and one always has 
the relation (|2TT >. In 6D, e.g., the radiation recoil force in 6D is 

fL = 1/ (<*V + |M«" - 4« 2 ^) . (48) 

and the Schott terms is 

/schott = (" M + y« 2 « M + y +4«V +4(aa)vA , (49) 

the sum of two being orthogonal to the 6-velocity. The Schott terms is again given 
by the finite part of the integrated bound momentum. 

In odd dimensions one always have tail terms. Split of the retarded field into the 
self- and the radiative parts is always possible, and the substitution into the equations 
of motion leads to divergent and finite terms. Divergent ones can again be absorbed 
introducing the counterterms. The split of the stress-tensor into the emitted and the 
bound parts is more subtle. The situation is obscured by the fact that though the free 
field is massless and thus propagates along the null cone, the full retarded potential 
fills the interior of the past light cone. Still one is able to obtain the general formula 
for radiation momentum which is no more associated with the retarded proper time 
on the world-line fl22l . 



4 Local method for curved space-time 

An approach initiated by DeWitt and Brehme and further applied to linearized grav- 
ity in [23 1 appeals to computation of the integral of the stress-tensor in the world- 
tube surrounding the world-line. This is similar to our calculation of the bound mo- 
mentum. However, in curved space-time the split of the stress tensor into the emitted 
and bound parts becomes problematic, so the complete analysis of the balance be- 
tween radiation, the kinetic momentum and the bound momentum is not available. 
Meanwhile, to compute the total reaction force one can use much a simpler calcu- 
lation substituting the retarded field directly into the equations of motion [ 17]. This 
approach was also formulated in higher even-dimensional space-time in the paper 
Ifl9ll which we follow here. 
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4.1 Hadamard expansion in any dimensions 



As in D = 4 iflOl . the curved space Green's functions for massless fields in other di- 
mensions can be constructed starting with the Hadamard solution. For simplicity we 
consider here the scalar case. The scalar Hadamard Green's function Gh(x,x?) is a 
solution of the homogeneous wave equation □ a G//(x,jc / ) =0, where □ = gjivV I V v . 
The procedure consists in expanding Gh(x,x?) in terms of the Synge world func- 
tion a(x,x'). For D = 4 the Hadamard expansion contains two terms singular in a, 
namely, a -1 and ln<7. In higher dimensions one has to add other singular terms, 
and by dimensionality it is easy to guess that each dimension introduces an addi- 
tional factor a -1 / 2 . Thus, the Hadamard expansion in D — 2d dimensions (d ^ 3/2 
is integer or half-integer) generically must read 



Gh(x,x') = 



(27t) a 



(50) 



where g n = g„(x,x'), v = v(x,x / ) are two-point functions. It can be shown, that in 
odd dimensions we actually have only odd powers of a -1 / 2 , and in even dimensions 
— only even powers, that is, an expansion in terms of inverse integer powers of <7. 
The logarithmic term is present only in even dimensions. 

Substituting d50b into the wave equation, in the leading singular order we will 
have: gD = A l l 2 . In the next to leading order we obtain the equation: 

2 dtfo- x o»+g D - X Uo-{D-\ )g D _ i = 0. (51) 

which does not have analytic solutions, so go-\ = 0. For D = 3 this means the ab- 
sence of the logarithmic term. Similarly, considering the equation for gm-i-2k): k e 
N we find go-i-ik = 0. This means that for an even-dimensional space-time the 
Hadamard Green's function contains only integer negative powers of a plus log- 
arithm and a regular part, while in the odd-dimensional case - only half-integer 
powers of a plus a regular part. 

For the sequence of Green's functions in the flat space-time, the one in D + 
2 dimension is proportional to the derivative of the Green's function in the twice 
preceding dimension D. In fact, in even dimensions the symmetric Green's function 
is the derivative of the order d — 2 of the delta-function: G D ~ S d ~ 2 (<r), a = 
(x — x') 2 /2 and thus, G D+2 ~ dG D /do. Applying regularization, S((x — x') 2 ) = 
lim e ^ + o5(|(x — x l ) 2 \ — e 2 ), we obtain 

G D+2 - dG D /de 2 . (52) 

This relation has a consequence that the Laurent expansion of the Lorentz-Dirac 
force in terms of e in the even-dimensional Minkowski space has only odd negative 
powers, and no even terms. So the number of divergent terms in the self-action in- 
creases by one for each next even dimension. In curved space passing to the standard 
notation notation for g n we have 
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^-^(|^I + vlna + wj, (53) 

where 

«o = ^ 1/2 

and we denoted v = w^-ii w = uj. Applying □ with respect to 
x^, we obtain the system of recurrent differential equations for Ui(x,x') . Integrat- 
ing them along the geodesic connecting the points x,xf, one can uniquely express 
U\{x,x!) through uq(x,x'). Furthermore, U2 is expressed through u\, etc. 

To relate the coefficient functions uf(x,x') in different dimensions we first ob- 
serve that mo = A 1 / 2 for any D. It is worth noting, however, that in the expansion in 
terms of a 

ug = 4 1/2 = l + l/12R ap a a a p + ... (54) 

where the tensor indices a, /3 run the values corresponding to D. With this in mind, 
we can write = for any D, D' . The next equation in the recurrence gives for 
D,D' > 5 

jjf _ D d-2 



L1 1 — Mi — ; . 

1 1 dl - 2 
Similarly for D, 7, D' = D + 2 one obtains: 



(55) 



"? +2 = "2^Y- (56) 



d-3 



Continuing this process further one finds 



d—\ da 

where the "direct" part of the Hadamard function is 

1 d-2 

1 U k 



1 Ar D 

d ° H = (2k)G^ d+ \ (57) 



G^dir-^pl^zr^. (58) 



4.2 Divergencies 

For the retarded Green's function one finds: 

G ret = ——— rT 0(x ,E (x) 2. i vfl(-ff) . (59) 

The first term constitutes the direct part of the retarded function with the support on 
the light cone: 

Gdir - 2(2^® [X] £ ml ' (60) 
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Similarly to the recurrent relations of the previous section, we obtain for the 
retarded Green's functions: 

d -^ = -2n{d-l)G°+ 2 . (61) 

Using the regularization 8(a) = lim £ ^ + o 8{o-S), where S = e 2 /2, we obtain for 
the direct part of the reaction force 

M dir 1 dfp dl \s,£) ff -r,s 

fo+2- D _ x ds ■ 

The limit S — > +0 has to be taken after the differentiation. The direct force is due 
to the light cone part of the retarded Green's function. This is not the full local 
contribution to the Lorentz-Dirac force. An additional contribution comes from the 
differentiation of the theta function in the tail term v9(o). In the scalar case this 
contribution vanishes, but in the electromagnetic case an extra local term arises: 

/loc = e " (fow]*v - [vva]z„) z v z a , (63) 

where the coincidence limit [v va ] depends on the dimension. The remaining contri- 
bution from the tail term will have the form of an integral along the past half of the 
particle world line. 

The direct part of the Lorentz-Dirac force contains divergences. To separate the 
divergent terms one can use the decomposition of the retarded potential suggested 
in the case of four dimensions by Detweiler and Whiting l24l l25l . In higher even 
dimensions we can follow essentially the same procedure. We define the "singular" 
part Gs of the retarded Green's function as the sum of the symmetric part (self) and 
the tail function v as follows 



v(x,x') . v(x,x') G(— a) 

Gs[X,X ) = G seK (x,X ) + 4 ( 2n y-i = G S elfdir(*,-* ) + ^^d-X ' ^ 



Here the direct part of the self function means its part without the tail v-term. The 
remaining part of the Green's function Gr(x,x!) = G^XjX 1 ) — Gs(x,xf) satisfies a 
free wave equation and is regular. Taking into account Qv = 0, it is clear that Gs 
satisfies the same inhomogeneous equation as G se if. The v-term in the second line 
of Eq. (|64| | is localized outside the light cone. Therefore the corresponding field (for 
instance, scalar), at an arbitrary point x will be given by 

T adv 

<h(x) = se lfdir+ 4( ^°y;-l J v( x ^(r))dT, (65) 



where the retarded and advanced proper time values T re t (x) , T a d v (x) are the intersec- 
tion points of the past and future light cones centered at x with the world line. 
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Inserting d65l l into the Lorentz-Dirac force defined on the world-line x = z(t), 
we observe that the integral contribution from the tail term vanishes and so the di- 
vergent part is entirely given by S eif dir - Using for delta-functions the point-splitting 
regularization we find that all divergent terms arise as negative powers of e. To prove 
the existence of the counterterms we consider the interaction term in the action sub- 
stituting the field as the S-part of the retarded solution to the wave equation. In the 
scalar case we will have: 

1 



5s = ^ J G s {x 1 x')p(x)p{x l )^g(z)g(z')dxdx l , (66) 

where a factor one half is introduced to avoid double counting when self-interaction 
is considered. Substituting the currents we get the integral over the world-line. Since 
the Green's function is localized on the light cone we expand the integrand in terms 
of the difference t = T — t' around the point z(t): 

5s ~ fdx [ Y,B kl (T)8M(t 2 -e 2 )t'dt. (67) 

J J k,l 

Here the coefficients By(z) depend on the curvature, while the delta-functions are 
flat: S( k > (t 2 — e 2 ). By virtue of parity, the integrals with odd / vanish, so only the odd 
inverse powers of e will be present in the expansion. Moreover, once we know the 
divergent terms in some dimension D, we can obtain by differentiation all divergent 
terms in D + 2, except for 1 /e term. The linearly divergent term corresponds to 
/ = 2k. The integral is equal to 







8® 



In four dimensions this term is unique. Applying our recurrence chain we obtain the 
inverse cubic divergence in six dimensions and calculate again the linearly divergent 
term. Thus in D = 2d dimensions we will get d— 1 divergent terms from which d — 2 
can be obtained by the differentiation of the previous-dimensional divergence, and 
the linearly divergent will be new. This linearly divergent self-action term in the 
action will have generically the form 

(-D i /•y (-i)*(2t-i)ii 

5 s ~ - / L BkMvdx. 

fc J k=o z 

Here the coefficient functions are obtained taking the coincidence limits of the two- 
point tensors involved in the expansion of the Hadamard solution. They actually 
depend on the derivatives of the world line embedding function z(t) as well as the 
curvature terms taken on the world line: 

B*,2*(t) - B k:2k {z,z, ...,R(z(t)),R^ v (z(t)), ...). 
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The vector case is technically the same, now one has to expand in in powers t in 
the integrand of 

S S = \J G^ ia (x,x')j >i (x)j a (Y ) s/g(zjg{^) dxdx'. 

From this analysis it follows that in any dimension the highest divergent term 
can be absorbed by the renormalization of the mass as in the generating four di- 
mensional case. To absorb the remaining d — 2 divergences one has to add to the 
initial action the sum of d — 2 counterterms depending on higher derivatives of the 
particle velocity. Typically the counterterms depend on the Riemann tensor of the 
background. 



4.3 Four dimensions 



In four dimensions the scalar retarded Green's function contains a single direct term 
localized on the light cone and a tail term: 



G iel = ^-6[£(x),x'} \A l l z d{o)-v6{o) 

47Z L 



(68) 



The retarded solution for the scalar field reads 



re tM = m q / -A l/2 8(a) + v9(a) 



dx' . 



(69) 



Differentiating this expression we obtain (j> v = d v (j) on the world line: 

T 

0v(z(t))=»k# J [-A xl2 8'{o)o v - Av /2 5(<t) -v8{a)a v +v v ]dr' , (70) 

CO 

where integration is performed along the past history of the particle. All the two- 
point functions are taken on the world-line at the points x = z(t) (observation point) 
and z' = z(t') (emission point). 

To compute local contributions from the terms proportional to delta-functions 
and its derivative, it is enough to expand the integrand in terms of the difference 
s = At = z — x 1 around the point z(t). The Taylor (covariant) expansion of the 
fundamental biscalar (7(z(t),z(t')) is given by l26l : 

CT(Z(T),Z(T')) = £ V(T(T, T)(T - x'f, (71) 
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where we denote as Q(z,t') the quantity Q(z(t),z(t')) for any Q taken on the 
world-line, and D is a covariant derivative along the world-line (also denoted by 
dot): 

Da = 6 = <T a z a , D 2 a = (J = G a pz a ^ + <J a z a , 

etc. Such an expansion exists since the difference s = X — %' is a two-point scalar 
itself: this is the integral from the scalar function / (— z 2 ) l ^ 2 ds, along the world-line 
from z(t) to z(t'). Taking the limits and using z 2 (t) = — 1, we find: 

<j(s) = -j-z 2 (T)^ + 0(s 5 )- (72) 

To obtain an expansion of the derivative of a overz^(T) one can expand C7^(t, T— s) 
in powers of s. This quantity transforms as a vector at z(t) and a scalar at z(t'), this 

a" (s) = * (V - z^ + z " ~) + 0(s% (73) 

where the index ji corresponds to the point z(t): = do(z,Z , )/dzu- Recall that 
the initial Greek indices correspond to z(f'). The expansion of 5(— cr) will read: 

8(-a) = S(s 2 /2) + s 4Z ^8'(s 2 /2) + ... (74) 

where the derivative of the delta-function is taken with respect to the full argument. 
Since the most singular term is A l l 2 8'{o)o v , the maximal order giving the non- 
zero result after the integration is s 3 . (Note, that in order to use our dimensional 
recurrent relations to obtain a reaction force in higher dimensions we should perform 
an expansion up to higher orders in s.) Thus, with the required accuracy, 8{—a) = 
28{s 2 ) 1 and all the integrals for the delta-derivatives are the same as in the flat space- 
time. This allows us to use the same regularization of the delta-functions with double 
roots 8{s 2 ) by the point-splitting. Expanding the biscalar A l l 2 and its gradient at z 
we have: 

4 l/2 = l+ S \^f^ dyA l/2 = S Rvxf (?5) 
12 o 

Combining all the contributions we obtain finally for the field strength on the world- 
line: 




After the renormalization of mass one recovers the familiar equation. 

Similarly, in the electromagnetic case one fins the retarded Maxwell tensor on 
the world line x = z(t): 
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F ret 
1 p.v 



= e I [uua; V 8(a)+Uua(T v 5'(a)+v V a;u. + Vua^v8(a)-{n <-> v}]z a dif, 

:(t) 



(77) 

Performing expansions on the world-line, one easily recover the DeWitt-Brehme- 
Hobbs result 



Jem — e 



2e 



(78) 

where n^ v = — z^z v - For geodesic motion the entire reaction force is given by 
the tail term. 



4.4 Self and rad forces in curved space-time 

Apart from the tail term, another new feature of the instantaneous momentum bal- 
ance between the radiating charge and radiation is the presence of the finite con- 
tribution in the self part of the reaction force originating from the half-sum of the 
retarded and advanced potentials. As we have seen in the previous section, in the 
case of Minkowski space the self part in four and other even dimensions is a pure 
divergence which has to be absorbed by renormalization of the mass and (in higher 
dimensions) the bare coupling constants in the higher derivative counterterms. In 
curved space, as was first shown in iflOl . the tail term in the equation for radiating 
charge moving along the geodesic with non-relativistic velocity (in weak gravita- 
tional field) contains apart from the dissipative term also the conservative force. 
This conservative force was later found for a static charge in the Schwarzschild 
metric fl4l . Here we would like to explore the significance of this result in view 
of the above analysis of self/rad decomposition. In the paper fT3l the splitting of 
the retarded field into the self and rad parts was not used. Considering the geodesic 
motion we have the equation: 

mz a =e 2 z p f v^V)^'- (79) 
Splitting the tail function with respect to T-parity 

C/}/ = V™!^/ + Ci/j/, (80) 

and repeating the calculations along the lines of lfl3l one finds in the weak-field 
non-relativistic case the corresponding (spatial) parts of the reaction force in terms 
of the flat space theory: 

fself = fdiv +fwWSM, frad = fsdiott; (81) 
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where the divergent term is the same as in the fiat space. The two finite terms 



are the WWSW force and the Schott force. Therefore, the WWSW force is a finite 
part of the T-even (self) contribution to the tail, while in the T-odd (rad) part repro- 
duces the Schott term which is precisely the same as in the flat space. Similar result 
holds for the scalar radiation. 



5 Gravitational radiation 

Gravitational radiation in the framework of linearized gravity on the curved back- 
ground may seem similar to scalar or electromagnetic radiation, but the similarity is 
incomplete. The difference is that radiation of non-gravitational nature emitted by 
bodies moving along geodesies in the fixed background is not influenced by the non- 
linear nature of the Einstein equations. For gravitational radiation this is not so. The 
scalar or vector linear field equations in curved space imply vanishing of the covari- 
ant divergence of the field stress tensor of matter field with respect to the background 
metric. In the gravitational case one has the Bianchi identity which generally does 
not imply the covariant conservation of the matter perturbation stress-tensor alone 
unless the background is vacuum. In most of the literature on gravitational radiation 
reaction the background is assumed to be vacuum. This seems to be satisfactory for 
the Schwarzschild or Kerr metrics. But physically we have to deal not with an eternal 
black hole, but with the collapsing body which is not globally vacuous. Meanwhile 
the conservation laws which may hold in the asymptotically flat case has to be con- 
sidered globally. It turns our that we must take into account contribution from the 
(perturbed) source of the background field as well. This is directly implied by the 
Bianchi identities. 



5.1 Bianchi identity 

B 

Let the background be generated by the stress-tensor T^ v . We are interested by grav- 
itational radiation emitted by a point particle of mass m moving along the geodesic 
of the background. Perturbations caused by the particle are assumed to be small so 
the particle stress-tensor 



fwWSM 



GMe 2 



fschott — r 




(82) 




(83) 
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B 

is the first order quantity with respect to T^ v ■ By construction, the tensor < f83b is 
divergence-free provided the particle follows the geodesic in the space-time. Since 
it is the first order quantity, it must be divergence-free with respect to the background 
covariant derivative up to the terms of the second order. Expanding the full metric 
§iiv = gp.v + f/jjiv, where K 2 = %%G, and the Einstein tensor 

£ M v = £ M v + (g4) 
we could naively expect the full Einstein equations to be of the form 

£ M v = K 2 ^v^j (g5) 

B B 

Since G^ v = K z T^ v , we then should have 

G^ = k 2 T^. (86) 

But the left hand side of this equation is not divergence-free with respect to the 
background covariant derivative. Expanding the full covariant derivative as 

V M =V M +V M (87) 

B B 

and taking into account the Bianchi identity for the background V ^ G^ v =0we 
obtain in the first order 

B 1 IB 

V M G MV = -V M G MV . (88) 

Therefore 

B m IB 

v M r MV = -v M r MV , (89) 



where the right hand side is the first order quantity. Thus Eq. ( I86I I is contradictory. 
Physically the reason is that we have to take into account the perturbation of the 
background 8T^ V caused by the particle, so that the correct equation should be 

G ^ v = k 2 (f Vv + 8T^ . (90) 

But this is not an equation for h^ v , since in order to find 8T^ V one has to consider 
the matter field equations for the background metric. The problem thus becomes 
essentially non-local. This non-locality does not reduce to that of the tail term in the 
DeWitt-Brehme equation. 



5.2 Vacuum background 

B 

If =0 the above obstacle is removed and one can proceed further with the 
linearized equations for h^ v . The derivation of the reaction force initiated in l27l 
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and completed in ||23l was based on the DeWitt-Brehme type calculation involving 
the integration of the field momentum over the small tube surrounding the particle 
world-line. As was noted later [28 1, this derivation had some drawbacks (for more 
recent discussion and further references see f29\ . One problem consisted in com- 
puting the contributions of "caps" at the ends of the chosen tube segment which 
were not rigorously calculated. Another problem was the singular integral over the 
internal boundary of the tube which was simply discarded. In addition, the usual 
mass-renormalization is not directly applicable in the gravitational case: due to the 
equivalence principle the mass does not enter into the geodesic equations. 

In ifTTl the local derivation of the gravitational reaction force was given which 
is free from the above problems and, in addition, is much simpler technically. It 
deal with the quantities defined only on the world line and does not involve the am- 
biguous volume integrals over the world-tube at all. The elimination of divergences 
amounts to the redefinition of the affine parameter on the world-line. 

We start with reparametrization invariant form of the particle action introducing 
the einbein e(x) on the world line acting as a Lagrange multiplier: 



S[z",e]=-l 



dx. (91) 



Variation with respect to and e{x) gives the equations 

and we obtain the geodesic equation in a manifestly reparametrization invariant 
form 



D 



(93) 



Assuming now that the particle motion with no account for radiation reaction is 
geodesic on the background metric, the perturbed equation in the leading order in K 
will read 



= \ [f V -^P) ( h x P -, v -2h vX ., p )z l f, (94) 

where contractions are with the background metric. 

The particle energy momentum tensor in our formulation will read 

8(x-z(x)) 



[e(xmr)iHx)^m d x, (95) 



and we choose the non-perturbed ein-bein eo =const as a bare parameter. After the 
calculation similar to that of the preceding section we obtain 
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Renormalization of the einbein is 



2" = 



(96) 



(97) 



Finally we choose the renormalized affine parameter so that z 2 = — 1 , which is equiv- 
alent to set e = m and obtain the MiSaTaQuWa equation. As was shown in [ 17 1 this 
equation remains valid in a class on non-vacuum metrics, in particular, for Einstein 
spaces. 



5.3 Gravitational radiation for non-geodesic motion 



^ = ~^b" v -ft , )i'» + tail. (98) 



If the particle world-line is non-geodesic, the radiation reaction force contains a 
putative antidamping term which is a local part of the rad contribution to the self- 
force: 

11k- 2 

T 

The reason is simply that the source of gravitational radiation is incomplete and 
the stress tensor is not divergence-free as required. Indeed, if the force is non- 
gravitational, one has to take into account the contribution of stresses of the field 
causing the body to accelerate. For instance, to describe gravitational radiation of 
an electron in the atom, one has to add the contribution form the Maxwell field 
stresses (spatial components, non-relativistic motion): 



where 



~ . . . . m ■ ■ at. . 

ny> J = -K 2 GT' J , r j =T' J +T' J , (99) 



£ m a zizi8 3 (X a ), fo = -fg (100) 



a=l,2 



4k (X 2 X 2 fl 2 



and X' a = x' —z' a (t). Using this source one can calculate the gravitational force and 
find the gravitational Schott term 

■ d 5 D'j mi m 2 

/Gshott = — TT-^r x b M = ~ — ■ < 101 ) 
Id at J mi +rri2 



where D !; is the quadrupole moment. 
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Note that this derivation of the Schott term is not based on the local calculation, 
the two-body treatment was necessary. This features seems to be general: gravita- 
tional radiation reaction form non-geodesically moving particle can not be described 
by some DeWitt-Brehme-type equation. 



6 Conclusions 

The purpose of this lecture was to discuss some subtle points associated with in- 
terpretation of the radiation reaction force. We have shown that the Lorentz-Dirac 
equation in classical electrodynamics describes the balance of three and not just two 
momenta: the mechanical momentum of the particle, the momentum of emitted ra- 
diation, and the momentum carried by the electromagnetic field bound to the charge. 
The total momentum is conserved, but this does not imply an instantaneous balance 
of the emitted momentum and that of the particle. The bound field momentum de- 
scribed by the Schott terms destroys the local balance. The total balance, however, 
is restored if one consider the situation when the charge has zero acceleration at 
the initial and final moments, or for a periodic motion subject to averaging. These 
considerations are equally applicable to radiation reaction of a charge in curved 
space-time and for gravitational radiation reaction. This explains, in particular, the 
necessity of averaging in calculating the evolution of the Carter constant in the Kerr 
field 02). 

A novel feature related to curved space is the existence of the finite WWSW 
force arising due to the tidal deformation of the bound electromagnetic field of the 
charge. This force is often interpreted as part of radiation reaction force, but one has 
clearly understand, however, that it has nothing to do either with radiation or with 
the Schott force. As we have shown, it is given by the T-even part of the retarded 
field, and thus present a finite remnant from the mass renormalization. 

Derivation of the reaction force of non-gravitational nature acting on a charge 
moving (both geodesically and non-geodesically) in curved space-time can be com- 
puted directly substituting the retarded field into the equations of motion, as is the 
Minkowski space. The regularization is easily achieved by the point-spliting, and 
divergences are eliminated by renormalization of mass. In higher dimensions one 
needs counterterms depending on higher derivatives of the velocity. Divergencies 
may contain the Riemann tensor of the background. 

Gravitational radiation reaction force can be obtained in a way similar to non- 
gravitational one only in vacuum space-time. In non- vacuum background the source 
of radiation apart from the local contribution from the particle must contain the 
contribution from the perturbed background. This can be seen from the analysis 
of the Bicnchi identity. This second contribution is non-local, so the possibility to 
obtain the equation of the DeWitt-Brehme type seems implausible. 

For a non-geodesically moving mass the formal derivation of the reaction force 
leads to putative antidamping effect. To cure this problem one has to take into ac- 
count the contribution of stresses forcing the mass to accelerate. Then in the non- 
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relativistic case one derives the gravitational quadrupole Schott term, but the deriva- 
tion is non-local. This is another example when the (quasi)local equation of motion 
with the reaction force does not exist. Here by quasilocality we mean the possibility 
of the tail term. 
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